
#pragma once

namespace cusp
{

template <typename LinearOperator,
          typename Matrix1,
          typename Matrix2>
void qr(LinearOperator& A,
        Matrix1& Q,
        Matrix2& R);

} // end namespace cusp

#include <cusp/detail/qr.inl>
